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In this paper, density functional theory calculations are used to explore the electronic and atomic 
reconstruction at interfaces between III-III/I-V oxides. In particular, at these interfaces, two di- 
mensional electron gases (2DEGs) with twice the interfacial charge densities of the prototypical 
LaTiOa/SrTiOa heterostructure are observed. Furthermore, a significant decrease in the band 
effective masses of the conduction electrons is shown, suggesting that possible enhancements in elec- 
tron mobilities may be achievable. These findings represent a framework for chemically modulating 
2DEGs, thereby providing a platform through which the underlying physics of electron confinement 
can be explored with implications for modern microelectronic devices. 
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Emergent phenomena at ABO3 oxide interfaces, 
e.g. two dimensional electron gases (2DEGs),[l] 
are paramount to understanding critical behav- 
ior arising from electron confinement; like metal- 
insulator transitions, [3 novel magnetic effects Q and 
superconductivity. 0, Q Recent theory and experiments 
have sought to identify the origin of 2DEGs at insu- 
lating oxide interfaces and surfaces in order to develop 
rules or concepts by which charge carrier densities and 
mobilities can be tuned. Growth under low O2 partial 
pressure has been shown to have (vacancy mediated) 
enhancements in conductivities that extend deep into 
bulk regions, making these materials undesirable for 
examining interfacial physics. @, 0] Conversely, annealing 
or growth under high O2 pressure eliminates oxygen 
vacancies resulting in consequential increases in resis- 
tivities; giving rise to intrinsic 2DEG behavior. 0, 0] 
In oxide heterostructures such as LaAlOs/SrTiOa and 
LaGaOa/SrTiOa this is understood in terms of polar- 
ization discontinuities at interfaces, i.e. the "the polar 
catastrophe" mechanism. 

@, [Ml Physically, the diver- 
gence of the electrostatic potential at atomically abrupt 
interfaces between insulating layered, charge-ordered 
III-III and a non-polar II-IV (usually SrTiOs) oxides is 
compensated for by charge accumulation (in this case 
an extra 1/2 electron per interface unit cell) at the 
interface. Interestingly enough, in these heterostructures 
this only occurs at TIO2 centered interfaces and requires 
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relatively thick layers of LaA103 or LaGa03. 

A chemically intuitive, ^-layer-doping, mechanism 
arises from the multivalent nature of transition metal 
cations, like Ti. For example, in LaTiOs/SrTiOs su- 
perlattices, the local environment of Ti cations next 
to a "dopant", LaO layer, splits the valence of Ti be- 
tween two possible charge states (+4 for SrTi03 and 
+3 for LaTiOs). Therefore, an equal mixture of Ti va- 
lence states (3+ or 4+) can be thought to reside at the 
interface giving an average valence of 3.5.[l4j-ll6j This 
has been confirmed through electron energy-loss spec- 
troscopy (EELS) measurements [T3| in which the dis- 



tribution of Ti 3+ cations away from the interface were 
found to be in good agreement with theoretical and ex- 
perimental carrier density profiles, fil Il8- 22 [ In a perfect 
system, the extra 1/2 electron, relative to Ti in SrTi03, 
defines the intrinsic limit of 2DEG carrier densities. [231] 
This electronic reconstruction is a hallmark of the ob- 
served two-dimensional conductivity and is accompanied 
by polar distortions [3, HH [25| (atomic displacements 
of the cations away from the interface) which effectively 
screen the electrons near the interfaces. Furthermore, 
the most commonly studied 2DEGs derived from II-IV 
(e.g. Sr 2+ Ti 4+ 03) oxides have interfacial carrier den- 
sities which are intrinsically limited to 0.5 e _ /interface 
unit cell. 

In this paper, density functional theory (DFT) calcu- 
lations are used to investigate the charge rearrangement 
at interfaces between I-V/III-III perovskites. The goal 
is to apply the above insights to control the interfacial 
charge density by increasing the intrinsic limit of elec- 
trons at an oxide heterointerface. The guiding principle 
is that the incorporation of a multivalent cation at an in- 
terface, where its desired valence states are +3 and +5, 
should allow for an average valence of +4, thus increas- 
ing the limit of extra interfacial charge to 1 e~ /interface 
unit cell. By examining heterostructures comprised of 
I-V/III-III oxides, it is confirmed that a total of 1 elec- 
tron per interface unit cell now populates the conduction 
bands (twice that of a corresponding LaTi03/SrTi03 su- 
perlattice). Similar to previous observations of 2DEGs, 
these electrons are primarily confined to t2 g orbitals on B 
cations near LaO interfaces, decay quickly into the bulk 
and are accompanied by large polar ionic distortions. In 
addition, calculated decreases in electron band effective 
masses suggest that improved carrier mobilities may be 
achievable. 

DFT calculations using the local density approxima- 
tion with a Hubbard U (LDA+U) 26} and ultrasoft 
pseudopotentials [27} as implemented in the Quantum 
Espresso simulation package [H, [29| were performed to 
study 1 LaXOs/7 KXO3 superlattices, where X=T& and 
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FIG. 1. DOS for (a) bulk SrTiOs, (b) 1 LaTi0 3 /7 SrTi0 3 , 
(c) 1 LaNb0 3 /7 KNb0 3 (d) 1 LaTa0 3 /7 KTa0 3 and (e) bulk 
KTa0 3 . All energies are relative to the Fermi level, Ep. (The 
scale for the heterostructures is 5x that of the bulk struc- 
tures.) 
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FIG. 2. Electronic band structure for (a) 1 LaTi0 3 /7 SrTi0 3 , 
(b) 1 LaNb0 3 /7 KNbOs and (c) 1 LaTa0 3 /7 KTaOs empha- 
sizing the partially occupied states near the Fermi surface. 



Nb. All superlattice calculations employed an 80 Ry cut- 
off and an 8 x 8 x 1 k-point mesh. Comparisons were made 
with a prototypical 2DEG system, 1 LaTi0 3 /7 SrTi0 3 . 
In all calculations, the in-plane lattice constants were 
constrained to the theoretical value of the majority com- 
ponent (i.e. KAO3 or SrTiOs) and the out-of-plane, 
c, lattice vector was optimized within the P4mm space 
group with lxl in-plane periodicity. Simultaneously, all 
ionic coordinates were relaxed until all Hellman-Feynman 
forces were less than 8 meV/A. The computed bulk 
KNb03, KTa03, and SrTiOs cubic lattice constants were 
3.951 A, 3.945 A, and 3.855 A, respectively. (Note: these 
values were obtained using standard LDA, i.e. without 
the inclusion of a Hubbard U). These are in typical LDA 
agreement with experimental values of 4.000 A, 3.988 A 
and 3.901 A, respectively. For all heterostructure cal- 
culations, a Hubbard U=5 eV for 5-cation d-states was 
found to be appropriate. Similar U values were used 
in previous studies of LaTiO 3 /SrTiO 3 .0 H3 Band ef- 
fective masses were compute using quadratic fits of the 
partially occupied bands. 

Figure Q] depicts the density of states (DOS) for bulk 
SrTi0 3 , bulk KTa0 3 , 1 LaTi0 3 /7 SrTi0 3 , 1 LaNb0 3 
/7 KNb0 3 and 1 LaTa0 3 /7 KTa0 3 superlattices. Both 
SrTi0 3 and KTa0 3 , using LDA (i.e. no Hubbard U) 
have relatively large band gaps of 1.7 and 1.8 eV, respec- 
tively. Although smaller than experiment, these are con- 
sistent with LDA's underestimation of oxide band gaps. 
In agreement with previous studies, 1 LaTi0 3 /7 SrTi0 3 
has occupied states just below the Fermi level, Ep , that 
sum to 1 electron (or rather 0.5 e~/interface unit cell). 
The electronic band structure plot (Fig. [2^) indicates 
that they directly contribute to transport (i.e. cross the 
Fermi level). An analysis of the orbital projected DOS 
indicates that these states are derived mainly from Ti ti g 



states, with the two lowest energy, light electron, bands 
coming almost entirely from d xy orbitals on the inter- 
facial Ti ions and the remaining occupied states being a 
mixture of t^g states on all of the Ti cations. More impor- 
tantly, in 1 LaNb0 3 /7 KNb0 3 and 1 LaTa0 3 /7 KTa0 3 
these occupied electronic states sum up to exactly 2 elec- 
trons (i.e. 1 e" / interface unit cell). An examination 
of the electronic band structure shows that they con- 
tribute to the Fermi surface (see Fig. [2)3 and c). Orbital 
projected DOS indicate that they are derived mainly 
from the 5-cation d-states (Nb/Ta) with dominant elec- 
tronic contributions arising from partially occupied d xy 
orbitals of 5-cations at the LaO interface. In all three 
heterostructures light electronic bands crossing Ep are 
parabolic around T and heavy bands extend along the 
T-X direction. This is a characteristic feature of 2DEGs 
and is consistent with recent angle-resolved photoemis- 
sion spectroscopy (ARPES) results for 2DEGs at SrTi0 3 
surfaces. 2^, 3(| 



Figure^ displays the spatial distribution of the con- 
duction electrons (i.e. arising from states between Ep 
and ~ -1.8 eV) as a function of distance away from 
the LaO layer. Similar to previous theoretical and ex- 
perimental results for 2DEGs at heterointerfaces and 
the SrTi0 3 surface we find a build-up of charge of 
roughly 0.22 electrons on interfacial Ti cations in the 
1 LaTi0 3 /7 SrTi0 3 superlattice. jEHH-El] (Note: the to- 
tal atom projected DOS adds up to 0.95 e~ and is scaled 
to unity in the plots.) This charge quickly decays to 0.06 
electrons in the center of the slab, indicating a decay 
length of roughly 3-4 unit cells. On the other hand, we 
observe that the LaNb0 3 /KNb0 3 and LaTa0 3 /KTa0 3 
heterostructures have a build-up of 0.52 and 0.60 e~s on 
the interfacial Nb and Ta ions, respectively, (again the 
projected DOS only sums to 1.8 e~s and is uniformly 
scaled to 2.0). Surprisingly, discernible differences in the 
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Relative z-coordinate 

FIG. 3. (color online) [top] Representative 2-dimensional pro- 
jection of a relaxed superlattices. (a) Charge distribution and 
(b) magnitude of B- and (c) j4-cation off-centering as a func- 
tion of relative 2-coordinate for the superlattices studied. All 
distances are relative to LaO planes. 



decay of charge away from the LaO interface are seen. 
In LaNb0 3 /KNb0 3 , there is still roughly 0.15 e~/unit 
cell area in the bulk region, whereas the charge density 
in LaTa0 3 /KTa0 3 drops off much more sharply, falling 
to less than 0.1 e~/unit cell area. This deviation may 
be linked to differences in dielectric constants. Theoret- 
ical calculations have demonstrated that a larger dielec- 
tric constant can result in a greater spread of electrons 

KTaOsand SrTi0 3 have 
" El while KNb0 3 



away from the interface. [311132 
very similar dielectric constants, |3S 
has a much stronger dependence of dielectric constant on 
phase. [35[ In fact, KNb0 3 can have a dielectric constants 
that is a few orders of magnitude greater than SrTi0 3 
and KTa0 3 , making the observed behavior reasonable. 
Regardless, these results clearly indicate a significant en- 
hancement in the concentration of electrons near the in- 
terface for III-III/II-IV heterostructure relative to the 
III-III/I-V system, with thelower dielectric constant ma- 
terials being more conducive to confining electrons to the 
interface. 

In accordance with the observed interfacial electronic 
reconstruction we find considerable polar distortions of 
the A- and 5-cations away from the LaO layer. These 
off-center displacements gradually return to zero in the 
center of the SrTi0 3 /KX0 3 layers. (See Fig. 1 [top] for 



TABLE I. Structural parameters and relative effective masses 
for the heterostructures studied. a , c/a and m e /m denote 
the cubic lattice parameter, c/a ratio and the relative effec- 
tive masses of the two lowest energy partially occupied bands, 
respectively. Note: the curvature of the electronic structure 
around F is symmetric, i.e. the effective mass along the F-X 
and T-M directions are essentially equal and thus only one 
value is reported for each band. 



System 


Oo [A] 


c/a 


m e 
1 


/m 
2 


1 LaTi0 3 /7 SrTiOs 


3.885 


8.117 


0.49 


0.59 


1 LaNb0 3 /7 KNb0 3 


3.951 


8.085 


0.35 


0.41 


1 LaTa0 3 /7 KTa0 3 


3.945 


8.009 


0.30 


0.35 



a representative 2 dimensional projection of the atomic 
structure of a 1/7 heterostructure). Similar to previ- 
ous DFT results we find that the magnitude of SrTi0 3 
off-centering is 0.18 A and 0.13 A for the A- and id- 
eations near the interface, respectively. [3, [25[ Figure [3Jd 
and c show that the magnitudes of the B- and A- 
cation off-center displacements in the LaX0 3 /KX0 3 
structures are all appreciably larger than those in 
LaTi0 3 /SrTi0 3 . These larger polar distortions are con- 
sistent with the need for larger interfacial polarizations to 
adequately screen the interfacial charges. Unexpectedly, 
the LaTa0 3 /KTa0 3 superlattice while having the larger 
interfacial charge exhibits smaller polar distortions than 
LaNb0 3 /KNb0 3 . Although the difference in the com- 
puted interfacial charge may be a consequence of scaling 
of the charge densities, it may also be a direct indica- 
tor of the magnitude of charge screening in these two 
materials arising from the differences in their respective 
dielectric constants - with the lower dielectric constant 
material being more strongly screened. A second expla- 
nation is that KNb0 3 , unlike KTa0 3 , has a polar ground 
state which may be more favorable to inducing polar dis- 
tortions. In which case, the polar nature of the KNb0 3 
structure may be more advantageous as recent studies 
have proposed the use of an electric field as a switch for 
controlling interface conductivity. 24 , 3(| 

Finally, Table U lists the relative band effective masses, 
m e /m, of the two lowest energy partially occupied bands 
(see Fig. [3]). Remarkably, the computed m e /m values for 
SrTi0 3 , are in excellent agreement with ARPES mea- 
surements per formed for SrTi0 3 surface 2DEGs (0.5 - 
0.6 m„/m) .[30] While the origin of the 2DEGs at SrTi0 3 
surfaces is almost certainly due to a different mecha- 
nism than the <5-doped structures, this result is suggestive 
of characteristic electronic structure features. It should 
be pointed out that La doped SrTi0 3 has considerably 
higher effective mass (m e /m > 1) and that the band 
effective masses predicted here neglects correlation ef- 
fects. Moreover, the 1 LaX0 3 /7 KA0 3 structures were 
found to have significant decreases in m e for these bands 
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(which dominate the electron density at the interfaces), 
implying possible increases in carrier mobilities in the III- 
III/I-V superlattices. Of course, semiconductor physics 
tells us that increases in charge densities (as observed in 
the 1 LaX0 3 /7 KX0 3 based heterostructures) are often 
accompanied by decreases in carrier mobilities. Typi- 
cally increased carrier densities arise from increases in 
dopant concentrations, which act as scattering centers 
thereby reducing the mean free path of the carriers. In 
the above materials the LaO layers could be considered 
dopants (i.e. <5-doped layers). Since the concentration 
of La remains constant at the interface there may be 
little change in the scattering lifetimes, r. Hence, devia- 
tions in electron mobilities, fi e , may be more affected by 
changes in m e (where [i e = er/m e ). Furthermore, recent 
dynamical-mean-field calculations of oxide heterostruc- 
ture 2DEGs suggest that correlation effects are least af- 
fected by dopant concentrations and that the band effec- 
tive masses have the most deviations. (37l| 

In summary, using first principles methods, it 
is demonstrated that superlattices comprising I-V 
(K 1+ [Nb/Ta] 5 +0 3 ) and III-III (La 3+ [Nb/Ta] 3 +0 3 ) per- 
ovskites have twice the interfacial charge densities of III- 
III/II-IV(SrTi0 3 /LaTi0 3 ) superlattices. Here, the flex- 
ibility of multivalent cations, like Nb and Ta, leads to 
an intrinsic limit of 1 e _ per interface unit cell, twice 
the limit of previously studied III-III/II-IV and I-V/III- 



III [2J, |38( superlattices. Also, changes in electron ef- 
fective masses (with no change in dopant levels) imply 
that further enhancements in mobilities may be achiev- 
able. As such, this research highlights a viable path for 
enhancing the properties of 2DEGs. In addition, devia- 
tions in polar distortions and charge redistribution em- 
phasize the need for a better understanding of the rela- 
tionship between dielectric constant and electronic and 
atomic reconstructions. In agreement with other theo- 
retical models, 31, 32[ these results suggest that lower di- 
electric constant materials (e.g. KTa0 3 ) should be more 
ideal for realizing truly two-dimensional electronic gases. 
Alternatively, polar oxides like KNb0 3 may be more suit- 
able for applications where it is necessary to switch the 
charge carrier concentrations. Naturally, synthetic lim- 
itations related to effectively reducing ions like Ta 5+ to 
Ta 3+ may exist. A feasible route may be to substitute 
interfacial Ta cations with a cation, like V, that is more 
easily reduced. Ultimately, these results present a chem- 
ically intuitive framework (in the absence of factors such 
as O vacancies) through which intrinsic carrier concen- 
trations, and perhaps even carrier mobilities, of oxide 
heterostructure 2DEGs can be tuned and may be use- 
ful in device engineering (39j and in controlling quantum 
phenomena due to electron confinement. 
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